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We investigate two-component ultracold fermionic atoms with repulsive interactions trapped in 
an optical lattice with a ladder structure. By applying the Bogoliubov-de Gennes equations to an 
effective t- J model in the strong correlation limit, we discuss how the spatially-modulated spin-singlet 
pairs with d-wave like symmetry are formed in the systems with trapping potentials. Furthermore, 
a close examination of the condensation energy as well as the local average of potential, kinetic and 
exchange energies by means of the variational Monte Carlo method elucidates that local particle 
correlations enhance the stability of the superfluid state via substantial energy gain due to singlet 
pairing in the high particle density region. 
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I. INTRODUCTION 



Ultracold atomic systems have attracted much inter- 
est since the successful realization of superfluidity in al- 
kali metal atoms [H, 0, 0, 0, B El- Among a number 
of interesting issues addressed in this field, fermionic 
atoms in an optical lattice @, formed by loading 
the ultracold atoms in a periodic potential, currently 
provide an extremely hot topic, where the superfluid 
Q and the Mott insulating states [l(| as well as the 
Fermi surface [k0 | have been observed recently. Exten- 
sive theoretical studies on attractive fermions in optical 
lattices have treated a wide variety of remarkable phe- 
nomena, such as the BEC-BCS crossover [l2[ 13 L 3, the 
superfluid-insulator transition [l6l. 17l. Il8l. 191. and 
the possibility of a supersolid state rgaTSSTEa 123, SI- 
Furthermore, Fuldc-Fcrrell-Larkin-Ovchinnikov (FFLO)- 
type superfluid states have been suggested in systems 
with trapping potentials [H, [H, S3, 13, IH HI • 



Fermionic systems with repulsive interactions, which 
have also been studied extensively, provide an ideal play- 
ground where we could realize various intriguing phe- 
nomena expected for hi ghly correlated electrons in con- 



densed matter physics [31 



32J, 133J, 13J]. A remarkable 
advantage in atomic gases over condensed matter sys- 
tems is that their experimental parameters are highly 
controllable. For example, the depth of the optical lat- 
tice potential is controlled by the light intensity, which 
leads to the precise alternation of atom tunneling be- 
tween adjacent sites. Furthermore, one can easily create 
low-dimensional systems by tuning the laser beam con- 
figuration with the use of interference effects. From a 
viewpoint of quantum simulators of correlated electron 
systems, it may be particularly interesting to study the 
superfluidity in two-dimensional (2D) or ladder-type op- 
tical lattices since they are directly related to unconven- 
tional superconductors observed in 2D high-T c cuprates 
and also some oxides with a ladder structure 



35,|3J,|3j 

38l | . Therefore, it is an important issue to search for 



an unconventional superfluid state in repulsive fermions 
loaded in optical lattices. In this connection, we note that 



there is a crucial difference between the optical lattices 
and the ordinary solid state systems, i.e. the presence of 
trapping potential in the former system. It is thus de- 
sirable to clarify the effect of trapping potential in order 
to discuss the superfluid states of repulsive fermions in 
optical lattices. 

In this paper, we study the stability of a d-wavc like 
superfluid state of ultracold fermionic atoms in an optical 
lattice with a ladder structure. This study is motivated 
by the recent successful realization of a superlattice with 
double-well structure, which should certainly stimulate 
a systematic study of superfluidity on ladder-type opti- 
cal lattices in the near future. As mentioned above, it 
is important to take into account the effects of repul- 
sive interactions in the presence of inhomogencity due to 
a trapping potential. By using the standard mean-field 
technique, we first describe the spatially-modulated spin- 
singlet pairs in a trapping potential. We further employ 
the variational Monte Carlo (VMC) method to take into 
account particle correlations more precisely [3{J H(| ■ It is 
clarified how the superfluid state with d-wave like sym- 
metry is stabilized in our trapped system with particular 
emphasis on the role played by local particle correlations. 

This paper is organized as follows. We introduce 
the model Hamiltonian in £}TTl and investigate how the 
spatially-modulated singlet pairing with <i-wavc like sym- 
metry is realized in mill In £|IV[ we use the VMC method 
to discuss the stability of the superfluid state in the pres- 
ence of strong correlations. A brief summary is given in 
the last section. 



II. MODEL HAMILTONIAN 

We study a superfluid state of ultracold fermionic 
atoms in an optical lattice. Here we consider a two- 
component fermionic system with Fermi-Fermi mixture 
such as 171 Yb and 173 Yb [4l|. Alternatively, we can deal 
with fermions with two accessible hyperfine levels such as 
\F,m F ) = |9/2, —9/2), |9/2, —7/2) for 40 K [ll|, where F 
and mp are the total atomic angular momentum and its 
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FIG. 1: (Color online) Schematic pictures of a ladder struc- 
ture [L x 2 sites) . Vb is the depth of the harmonic potential in 
the system. Ai^+i and Ai,i+£ are the nearest neighbour pair 
potentials along the leg and the rung directions (see text). 



magnetic quantum number, respectively. By specifying 
the two different fermions by pseudo-spin indices a =|, J,, 
we can describe the ultracold atoms in an optical lattice 
in terms of the ordinary Hubbard Hamiltonian with on- 
site interactions [42], Hjj]. Here we focus on the strong 
coupling regime, where doubly occupied states hardly 
appear in the ground-state configuration, and then the 
antiferromagnetic superexchange terms between two ad- 
jacent sites become important to determine the nature of 
the ground state. In this regime, an effective model, the 
so-called t — J model, is more convenient to discuss the 
ground state properties. The model Hamiltonian thus 
reads, 

7~it-j = Tio + H , (1) 
W " = "*E (clc jtr + h.c) +^(^-M)n iCT (2) 



(3) 



where cj (c^) is a creation (annihilation) operator of 



a fermion with spin a at site i, 



and rii = 



riif + riii. S t = (1/2) J2 a ,p Cia&apCip, where cr is the 
Pauli matrix, t and </(> 0) are the tunneling matrix 
and the antiferromagnetic exchange interaction between 
nearest neighbour sites, and fj. the chemical potential. 
The interacting fermions are trapped in the harmonic 
potential Vj. Note that the doubly occupied states are 
excluded at each site in the Hamiltonian (fTJ) due to the 
strong onsite repulsive interactions. 

In this paper, we deal with the optical lattice with a 
ladder structure schematically shown in Fig[TJ Owing 
to the existence of the rung in the ladder, this quasi- 
1D system may have essence of 2D systems: namely a 
d-wave superfluid state may be realized in the strong 
correlation limit. It is important to note that this 
sort of ladder system can be experimentally realized 
by double well potentials [H|, |45[. The tunneling ma- 
trix t and the exchange coupling J are given such that 
t ~ (4/V^)Er(E r /v ) 3/4 cxp[-2v/V^] a nd J ~ 

within the harmonic approximation [34], [43 . [43| , where 
a s is the s-wave scattering length and E R = ti 2 k 2 /2m 
the atomic recoil energy in the optical lattice poten- 
tial, to, vq and k are the mass of the atom, the in- 



tensity and the wave number of the lattice potential 
Viati'x) = vo J2i sin(fcrf). Thus the exchange interaction 
J ~ t 2 /U can be experimentally controlled. Further- 
more, some elegant techniques to control J have been 
proposed by means of the spin-dependent lattice poten- 
tials [itl |47| . [48| and the superlatticc potentials [49| . 
Therefore, the ladder system has some advantages in 
studying the stability of the superfluid state with c?-wave 
like symmetry. In this study, we consider the ladder sys- 
tem with L x 2(L = 50) sites and fix the total number of 
atoms as = N± = 25, where N a = n ia- The har- 
monic trapping potential is given as V% = at four sites 
around the center of the system and Vi = 4. Of at the 
edges. Note also that V t = Vl-\-% = V l+L = V 2 L-i-i- 
We fix the chemical potential fx at the Fermi level in the 
non-interacting case and use t as an energy unit. 

In the following, we discuss the possibility of a e?-wavc 
superfluid state in our optical ladder system. Since the 
fermionic atoms arc confined in the inhomogencous po- 
tential, it is not straightforward to deal with the effects of 
particle correlations. In this paper, we make use of the 
Bogoliubov-de Gennes (BdG) equations and the VMC 
method with a proper trial function to deal with this 
problem. 



III. SPATIALLY-MODULATED SINGLET 
PAIRING STATE 

In this section, we discuss how the spatially-modulated 
superfluid state is realized by means of the BdG mean- 
field theory. We here introduce site-dependent singlet 
pairing mean-fields Ay = (c^Cjf) (= (cj|Ci|)). In terms 
of the Bogoliubov transformation Cj CT = ~^2 X {ufa\ a — 

av^a^g}, we obtain the following BdG equations [5Cj| . 



H F ■ 



with 



Hij = -*5<y) + (Vi - (j)S, 



I J ■■ 



F, 



and the self-consistent equations, 



A 



(4) 



(5) 
(6) 



(7) 



where 8uj\ is the Kronecker delta for the nearest neigh- 
bour sites i and j. We ignore the constraint for dou- 
bly occupied states in the Hamiltonian |T]) for a while, 
which largely reduces our numerical efforts. In the peri- 
odic case, it is known that the constraint mainly yields 
the renormalization of tunneling matrix t [50j . In our 
inhomogencous case, however, the situation is more com- 
plicated. The constraint will be carefully treated in the 
next section. 
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FIG. 2: (Color online) (a) Spatial distributions of the parti- 
cle density (rii)(= (th+l)), (b)-(f) the pair potentials along 
the leg direction Aj^+i (open circles) and the rung direction 
Ai, i+ L (filled circles) for J = 0.5,0.8,1.0,1.2 and 1.4 (see 
FigEE). 



We first calculate the site-dependent particle density, 
(jn) = ^2 X \vi\ 2 . The results are shown in Fig[2] (a). It 
is seen that the particles are smoothly distributed in the 
lattice. Its profile is mainly determined by the curvature 
of the harmonic potential, so that it is hardly affected 
by the exchange interaction J. On the other hand, the 
pair potential between the nearest neighbour sites A^- 
strongly depends on the magnitude of the exchange in- 
teraction, as seen in Figs[2] (b)-(f). For J < 0.5, the 
pair potential is too small to detect within our numeri- 
cal accuracy, where a singlet pairing state is not realized, 
and the system is still in the normal ground state. On 
the other hand, for J > 0.5, the pair potential takes fi- 
nite values, as shown in Figs[2] (b)-(f), implying that the 
intcrsite exchange coupling stabilizes the singlet-pairing 
supcrfluid state. Here, we find two characteristic fea- 
tures in the profile of the pair potential. One is the cusp 
(or peak) structures around i = 5 and i = 15, and the 
other is the plateau around the center. The above char- 
acteristic properties reflect atomic states inherent in our 
ladder system. To see this clearly, we introduce the in- 
tegrated local density of states (LDOS), W(i, J), around 



the Fermi level tF- The quantity is defined as, 

w(i,j) = [ F P (i,0<%, 

= e r + j wi^-go^' ( § ) 

where p(i, £) denotes the LDOS at site i, and £„ and 
4> n = {0o , <j>i, • ■ ■ 02L-i} are the eigenenergies and the 
eigenfunctions for the non-interacting Hamiltonian 7io, 
respectively. Let us first take a look at the LDOS p(i,£,) 
for Ho shown in Figlllll (a). It is seen that the shape of 
the LDOS at site i reflects the DOS for a uniform lad- 
der system that has a characteristic four-peak structure 
originating from the van-Hove edge singularity in one di- 
mension. At first glance, it is not evident whether this 
LDOS is really related to the spatially modulated pair 
potential. However, the integrated LDOS over the range 
of J, W(i, J), makes this point clear. In Fig lllll (b). it is 
found that the cusp (or peak) structures indeed appear 
around i = 5 and 15 in W(i, J), and the plateau with 
some fluctuations is formed around the center. This im- 
plies that the spatial variation of the pair potentials re- 
flect the integrated LDOS, in agreement with the fact 
that the particles around the Fermi level are relevant 
to the formation of spin-singlet pairs in the superfluid 
state. An important point to be noticed here is that 
the pair potentials along the leg and the rung directions 
have opposite signs as shown in Figs|2] (b)-(f), implying 
that the pairing state with c? a .2_ J) 2-likc symmetry is sta- 
bilized in our inhomogeneous system. This is naturally 
expected from the d-wave superconducting state realized 
in the doped ladder system without the trapping poten- 
tial [HI, [si, m, HI] . 

We have so far neglected the local constraint due to 
strong onsite repulsion at each site. Although this kind 
of approximation has been employed in some cases [50l . 
HI, Hy, [5?| , it is not sufficient to discuss the stability of 
superfluidity in our correlated system: for example, the 
particles in the normal phase arc completely free, and 
also the onsite pair potential An can be finite in the 
large J case, in contrast to the ordinary d-wave state in 
the homogeneous system. To remedy these pathological 
points, a careful treatment of local particle correlations 
is necessary, which will be done in the next section. 

IV. VMC STUDY OF PARTICLE 
CORRELATIONS 

In this section, we perform the VMC calculation to 
elucidate how the local particle correlations affect the 
supcrfluid state. Let us introduce a trial state, 

m=V G \$Bcs({&i j },g)), (9) 

where Vg = FIJI ~ n i\ n i\\ is a Gutzwillcr projection 
operator [58|], which completely excludes the doubly oc- 
cupied states at each site. \§BCs({Aij}, <?)) is a BCS 
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FIG. 3: (Color online) (a) LDOS p(i,0 = £ n \<t>?\ 2 5($ - £„) 
at i — 10 (dotted line) and i = 20 (solid line) :delta functions 
are broadened with width of r\ = 0.2t. (b) The integrated 
LDOS around the Fermi level W(i,J) in eq.® for J = 0.5 
(filled circles) , 0.8 (open circles) , 1.0 (filled squares), 1.2 (open 
square) and 1.4 (filled triangles). 



where g is a variational parameter controlling the ampli- 
tude of the pair potential, while A^- determines its spatial 
variation according to the solution of BdG equations. By 
optimizing the variational parameters so as to minimize 
the total energy, we discuss the ground state properties in 
the optical lattice system. It is still not easy to determine 
these parameters since the number of parameters is very 
large in the presence of the trapping potential. Recall 
here that in the solutions obtained by the BdG equa- 
tions, the spatial dependence of Ay is controlled by J. 
This observation enables us to approximately reduce the 
parameter space by regarding J in eq. as a variational 
parameter to control the distribution of {Ay} obtained 
from the BdG equations. This new variational parameter 
is denoted as J in the following. In order to incorporate 
the strong correlation effects in the order parameter, we 
here set the onsite component of the pair potential to be 
zero, Ajj = 0, and retain only the nearest neighbor com- 
ponent of pair potential, A y — > A^S^, which gives the 
most significant contribution to the superfluid state. This 
kind of approximation should be justified in the strong 
coupling regime. Note that the correlation effects in the 
normal state can be properly taken into account by the 
ordinary VMC procedure. 
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state with variational parameters ({A.y},g) to describe 
the superfluid state with spatially-modulated pairing cor- 
relations. This state is explicitly given as, 



I* 



BCS) 



oc 



U^n+Vnci^) |0), 

n 



|o>, 



with 



(10) 

(11) 
(12) 



where c ncr = Y^ifficia an< i N = 2Vj + iVj. The gap 
function A n is then given as, 



A n = -J(c nl C„ T ), 



(13) 



(c) J=0.3 J=0.5 
3=0.8 





FIG. 4: (Color online) Condensation energies as a function of 
the variational parameter g for (a) J = 0, (b)0.2, (c)0.3 and 
(d)0.4 obtained by the VMC method. The sets of {A 4j } used 
in each calculation are obtained for J = 0.5 (open squares), 
0.8 (filled squares), 1.0 (open circles), 1.2 (filled circles) and 
1.4 (open diamonds), which are shown in Figs[2] (a)-(d). 



By performing the VMC simulations with samplings 
~ 10 6 , we obtain the condensation energy AE = E super — 
Enormai, for the variational parameters (J,g) as shown 
in FigJU where E super and E norma i are the energies of 
the superfluid state with A„ ^ and the normal state 
with A„ = for all n. At J = 0, the model is reduced 
to the system only with correlated hopping. In the case, 
the formation of the singlet pairs gives rise to the loss of 
the kinetic energy. Therefore, the normal state is more 
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stable than the superfluid state with any values of J, as 
shown in FigfJ] (a). On the other hand, the increase in 
the exchange coupling J favours the superfluid state with 
singlet pairing correlations, where the normal state is no 
longer the ground state. In fact, it is found in Figs[4](b)- 
(d) that the condensation energy AE = E super — E n 

ormal 

has a minimum at a finite g although statistical errors 
are not so small in our calculations. This suggests that 
the spatially modulated superfluid state with d-wavc like 
symmetry is indeed stabilized by exchange interactions. 
It should be noted that the critical value of J ~ 0.2 is 
smaller than J ~ 0.5 estimated in the previous section 
without electron correlations, which implies that the cor- 
relation effects enhance the stability of superfluid state. 
At first glance, this conclusion seems a little bit unusual, 
since the correlation effects normally have a tendency to 
suppress the ordering. 

In order to elucidate the role of the onsite correlation 
effects from a microscopic point of view, we now compare 
the VMC results obtained with the Gutzwiller-projected 
state © with those of |*') = |$scs({Ay }, g)). Since 
the Gutzwiller projection is neglected in the latter, strong 
onsite correlation effects are not incorporated in the 
results obtained with \^f'). Here we compute the lo- 
cal average of potential, kinetic and exchange energies, 
which are respectively defined by Ef ot = Vi(rii), Ef™ = 

-*£*, (clc ja + h.c), Er = (J/2) Ej (Si ■ S 3 ->. 
Shown in Fig[S]are the profiles of particle density and the 
local average of energies obtained by the VMC simula- 
tions. It is seen that the onsite correlation effects hardly 
change the density profile of particles (FigO (a)). As 
a result, there is little difference in the potential energy 
Ef ot obtained with |*) and |*'), as seen in Fig[5] (b). On 
the other hand, the introduction of the onsite correlation 



effects significantly reduces E\ m around the center of the 
system (FigJS](a)), because the strong repulsive interac- 
tions renormalize the motion of particles, forming heavy 
fermion states in the region with high particle density. 
In fact, the ratio of the two kinetic energies evaluated 
by Gutzwiller-projected and non-projected states can be 
regarded as the site-dependent wave-function rcnormal- 
ization factor z%: ~ 1.0 (zj ~ 0.2) around the edge 
(center) of the system in Fig[5] (c). We note that these 
results hardly depend on J. On the other hand, it is seen 
in Figj5] (d) that the antiferromagnetic spin correlations 
are developed among the heavy particles around the cen- 
ter of the system, where the energy gain due to Ef x is 
enhanced in the presence of the strong correlation effects. 
Summarizing the above results, the onsite correlation ef- 
fects enhance the energy gain by the exchange interac- 
tions, although they cause the loss of the kinetic energy. 
An important point is that in the high density region 
((rii) ~ 1.0), the effective exchange coupling J divided 
by the renormalized hopping Z{t can be very large, which 
is efficient to stabilize the d-wavc superfluid mediated by 
antiferromagnetic correlations. This naturally explains 
why the critical value J ~ 0.2 obtained for correlated 
systems can be smaller than J ~ 0.5 for non-interacting 
systems. 

Note that the above conclusion holds true for the den- 
sity regime with n, < 1.0. For the system with higher 
particle density, the spatial distribution of particles is ex- 
pected to form a plateau around the center so as not to 
exceed n» = 1.0, because the strong particle correlations 
have a tendency to form the Mott insulating state. Even 
in such cases, the present conclusion may be applied to 
the strongly correlated conducting region in the vicinity 
of the Mott insulating state. 
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FIG. 5: (Color online) (a) Particle density {ru), and local 
average of energies: (b) potential energy Ef ot , (c) kinetic en- 
ergy E\ ln and (d) exchange energy Ef x obtained with the 
Gutzwiller-projected state |$) (open symbols) and the non- 
projected state \*$>') (filled symbols). We note that {ru), Ef ot 
and E\ m hardly depend on J. 



V. SUMMARY 

We have investigated the t-J model to discuss the sta- 
bility of the spatially-modulated superfluid state in a 
fermionic optical ladder system with harmonic confine- 
ment. By means of the BdG equations, we have clar- 
ified that the spatially-modulated superfluid state with 
d x 2_ y 2-\ihe symmetry can be realized via antiferromag- 
netic exchange interactions. It has been found that the 
spatial distribution of the pair potential has a character- 
istic profile with cusps and plateau, which reflects the 
local density of states integrated around the Fermi level. 
Furthermore, the VMC calculation of the condensation 
energy as well as the local average of potential, kinetic 
and exchange energies has shown that the local particle 
correlations enhance the stability of the superfluid state 
through energy gain due to singlet pairing in the region 
of high density filling ((n,) ~ 1.0). In particular, we have 
demonstrated that the strong renormalization of kinetic 
energy around the high density region plays an important 
role in stabilizing the singlet-pairing superfluid state. We 
think that this kind of behavior may be characteristic of 
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the systems with nonuniform distributions of the particle 
density. 

In this paper, we have exploited the t-J model as an 
effective model in the strong coupling regime. In this 
model the exchange interaction J is assumed to be free 
parameters, in contrast to the strong-coupling limit of the 
original Hubbard model. It thus remains an important 
problem to discuss the superfluidity by directly dealing 
with the Hubbard ladder system with a trapping poten- 
tial, which is now in progress. 

When the density of particles is increased, other states 
such as the Mott insulating state and the antiferromag- 
netically ordered state may appear around the center of 
the system (U [H, EHf EH . Therefore, it is particu- 
larly worth studying how the d-wave like supernuid state 
coexists or competes with the magnetically ordered state 
in the optical lattice system, which is to be addressed in 



the future study. 
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